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Abstract: Presented is an analysis of general scaling perturbations in a transmitting fiber. For 
elliptical perturbations, under some conditions an intermodc dispersion parameter characterizing 
modal PMD is shown to be directly proportional to the mode dispersion. 
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In the following paper we derive a generalized Hcrmitian Hamiltonian approach for the treatment of Maxwell 
equations in waveguides as well as develop a perturbation theory for the general class of scaling perturbations 
that include ellipticity and a uniform scaling of an arbitrary index profile. Because of the Hermitian nature 
of the formulation most of the results from the well developed perturbation theory of quantum mechanical 
systems can be directly related to the light propagation in the waveguides. Such formulation provides for 
an intuitive way of understanding PMD and birefringence in the elliptically perturbed fiber profiles. Region 
of validity of our theory extends to the case of large variations of the dielectric constant across the fiber 
crosscction and is limited only by an amount of re-scaling. Finally, we establish that if in some range 
of frequencies a particular mode behaves like a mode of pure polarization TE,TM (where polarization is 
judged by the relative amounts of the electric and magnetic longitudinal energies in a modal crossection) its 
inter- mode dispersion parameter r = | ^ | is related to its dispersion D as r = A(5|Z?|, where S is a measure 



of the fiber ellipticity and A/3 e is a split in a wavevector of a linearly polarized doubly degenerate mode of 
interest due to an elliptical perturbation. 

While there has been a wide amount of work done on estimating such quantities as local birefringence 
induced by perturbations in the fiber profile most of the treatments were geared toward understanding the 
low contrast, weakly guiding systems such as ubiquitous silica waveguides and are not directly applicable 
to the high contrast systems such as Bragg fibers, photonic crystal fibers and integrated optics waveguides 
which are steadily emerging as an integral part of the state of the art transmission systems. 

In deriving Hamiltonian formulation for the eigen fields of a generic waveguide exhibiting translational 
symmetry in longitudinal z direction we start with a well known set of Maxwell equations written in terms 
of transverse and longitudinal fields [1]. Assuming the form of the field 
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and introducing the transverse and longitudinal components of the fields as F = F t + F z , F z 
F t — (z x F) x z Maxwell equations can be rewritten as a generalized Hermitian eigen problem [2] 
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In this form operators on the left and on the right are Hermitian thus defining a generalized Hermitian eigen 
problem and allowing for all the convenient properties pertaining to such a form, including real eigenvalues 
(3 as well as orthogonality of the modes corresponding to the different /3's (for more discussion see [2]). 
Defining hermitian operator on the left of (2) as B and the one on the right A and introducing Dirac 



notation \ip >= 



E, 



we rewrite a generalized eigen problem as 
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with a condition of orthogonality between modes of (3 and (3' in the form 



< ipp\B\ipp, >= 8p,p>. 
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In the following we analyze uniform along z axis perturbations. If a perturbation such as general re-scaling 
of coordinates x sca i e( i = x(l + 5 x ),y sca i e d = y(l + 5 y ) is introduced into the system it will modify an operator 
A. A particular case of general re-scaling when S x = 5 y correspond to the uniform scaling of a structure, 
while the case of S x = —8 y corresponds to a uniform ellipticity. Denoting a correction to an original operator 
on the left of (2) SA, the new eigen values of the split doubly degenerate eigen mode are found by solving 
a secular equation [3] and gives 
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The inter-mode dispersion parameter being proportional to PMD [4] is defined to be equal to the mismatch 
of the inverse group velocities of the split due to the perturbation modes r = -4p \r which is, in turn, can 



be expressed in terms of the frequency derivative r 
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Now we derive a form of the perturbation operator for the cases of uniform scaling and uniform ellipticity. 
We start with an elliptical waveguide and a generalized Hermitian formulation (2) where the derivatives in 
operator A should be understood as the derivatives over the coordinates x sca i e d and y scaled- We then transform 
into the coordinate system in which an elliptical waveguide becomes cylindrical. Assuming normalization (4) 
after some combcrsome algebra [5] we arrive at the following expressions. 

Case of uniform scaling S x = 5 y = 5 



A/3 S =< ip0, m \8A\ipp tfn >= 26 j s ds 
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Another important result about change in the propagation constant of a mode under uniform scaling is 
that its frequency derivative is proportional to the dispersion of a mode. To derive this result we consider a 
dispersion relation for some mode of a waveguide ft = f(u>). From the form of Maxwell equations it is clear 
that if we uniformly re-scale all the transverse dimensions in a system by a factor of (1 + 5) then the new 
p = (3+ A/3 S for the same u> will satisfy (3 — ^^Cs • Decomposing the last expression in Taylor scries and 



collecting terms of the same order in S we derive expressions for A/3 S and its derivative A(3 S = S(uj^ 
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-X6D(u>) where D(u>) is a dispersion of the mode. 



Case of uniform ellipticity 6 X = —8 y = 5. A first order correction to the split in the values of propagation 
constants of the modes (/?, 1) and (f3, —1) due to the uniform re-scaling becomes (5) 
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A[3 e = 2| < V/3,iMlV73,-i > I = 2H J s ds 

2tjw| J s ds[(-e\E z \ 2 + \H Z \ 2 ) + 2Im(eE*E e - H*H e ) 
where E's and H's are those of the (/3, 1) mode. 

From expression (7) we find that the split between the degenerate modes due to the ellipticity is proportional 
to the difference in the longitudinal magnetic and electric energies in the crossection of a fiber. The rest of 
the crossterms in expression (7) usually do not contribute substantially to the split, unless special structures 
are considered where longitudinal magnetic and longitudinal electric energies are of the same order. 
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Fig. 1. Effective birefringence due to the elliptical perturbation in a double core high dielectric contrast 
fiber. Data is presented for the fundamental linearly polarized doubly degenerate mode. While the width of 
a ring was kept at a constant l/im corresponding inner radii R was varied in the interval (3,10)/xm. Split 
in a wavevector of an originally degerate mode due to the uniform elliptical perturbation of magnitude 8 as 
predicted by the perturbation theory (circles) is compared to the results of the Finite Difference numerical 
simulations (crosses). Excellent correspondence over the whole range of inner radii is observed. 

An important conclusion about PMD of a structure can be drawn when electric or magnetic longitudinal 
energy dominates substantially over the other (for a longer discussion see [5]). In the case of pure-like TE 
(E z ~ 0) or TM (H z ~ 0) modes split due to the uniform scaling (6) becomes almost identical to the split in 
the degeneracy of the modes due to the uniform ellipticity perturbation (7). Thus, in the case when the mode 
is predominantly TE or TM as judged by the amounts of the corresponding longitudinal energies in the 
crossection we expect A(3 S = A/3 e . As PMD is proportional to the intermode dispersion parameter r = 
and taking into account expressions for the frequency derivatives of A/3 S we arrive to the conclusion that for 
such a mode PMD is proportional to the despersion of a mode 
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We conclude by presenting the results of calculations of the normalized birefringence due to the uniform 
elliptical perturbation in the case of a double core high dielectric contrast fiber Fig.l. Fundamental doubly 
degenerate mode of m = 1 was studied in a ring-like fiber of 1.0 to 1.5 index contrast. Excellent correspon- 
dence between the predictions of the pertubation theory and, in principle, exact Finite Difference numerical 
simulations is observed. 
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